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Dust discs around intermediate mass and Sun-like stars in 
the 16 Myr old NGC 1960 open cluster 
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ABSTRACT 

We present an analysis of Spitzer IRAC (3.6-8/j,m) and MIPS (24/zm) imaging of 
members of the 16tj5°Myr old open cluster NGC 1960 (M36). Models of terrestrial 
planet formation indicate that rocky planets are likely to achieve their final masses at 
around 10-30Myr, and thus this cluster is at an interesting epoch for planet formation. 
We find 21 B-F5 type stars and 14 F6-K9 type stars which have 24/xm excess emission, 
and thus determine that > 30% of B-F5 type stars and > 23% of F6-K9 type stars 
in this cluster have 24/xm excess emission. These excess frequencies are similar to 
those observed in other clusters of similar age. Three early type stars have excesses at 
near-infrared wavelengths. Analysis of their SEDs confirms that these are true debris 
discs and not remnant primordial or transitional discs. None of the 61 sun-like stars 
have confirmed near-infrared excess, and we can place a limit on the frequency of 
8/xm excess emission around sun-like stars of <7%. All of the detected excesses are 
consistent with emission from debris discs and are not primordial. 

Key words: circumstellar matter - infrared: stars. 



1 INTRODUCTION 

Most young (1-2 Myr old) stars are surrounded by a gas-rich 
primordial disc with dust producing h igh levels of emission 
in th e infrared {Fir/F^ ^ 0.1, e.g. iKenvon fc Hartmannl 
Il995h . These are the birthplaces of planets. Current mod- 
els of extra-solar planet formation propose that dusty discs 
around a new star settle and km-s ized planetesimals aggre- 
gate on a short (<lMyr) timescale (jWeidenschilling &: Cuzza 
1993). The largest planetestimals undergo runaway accre- 
tion followed by oligarchic growth resulting in tens or hun- 
dreds of lO OOkm-sized b odies in their own cleared "feed- 
ing zones" |Klahrl [2008 ). These phases may take up to 
a few million years. Circumstellar gas accretes onto the 
star or onto large protoplanet ary cores to create gas giant 
planets (|Hartmann et alii 19981 ). Finally these planetary em- 
bryos collide and merge in a chaotic growth phase to form 
a few stable terrest rial planets over 10-100Myr (see e.g. 
IWeidenschillin3ll977l ). These or other pr ocess remove the 
primordial discs on tim e scales of 3-7Myr ijHernandez et al.l 
120071 : lHillenbrandll2008l : ICurrie et al.ll2009l h By 10 Myr al- 
mo st all remain ing discs are optically thin debris discs (see 
e.g. IWvatdl2008l and references therein). Emission from such 
discs arises from second generation dust populations pro- 
duced in collisions between planetesimals. These dust grains 
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absorb and re-radiate star light at wavelengths typically 
>10/xm. 

The Spitzer Space telescope has been used to study 
the evolution of debr is discs in a sta t istical manner (see 



e.g. iRieke et all 120051, ISu et all [20061, ISiegler et all |2007| . 
iRebull et all I2OO8L ICarpenter et al.l 120091 ). These studies 
have tried to answer the question of why two apparently 
similar stars can have very different levels of excess emis- 
sion. To date, the clearest dependency is on age. In A and 
early F-type stars there is evidence for a peak in the upper 
envelope of excess emission at 10-20Myr be fore a decay in 
proportion with time (see e.g. I Wvatd [20081 and references 
therein). For solar- type stars the number of observed ob- 
jects is smaller and so correlations are harder to establish. 
Based on current evidence the decay of 24/im excess around 
solar-type stars appears to follow a similar pattern to the 
A stars but on a timescale that is an order of magnitude 
shorter (a drop from 40% to 20% of stars with 24/im ex- 
cess occurs between 10 and lOOMyr for solar-type stars, and 
between 100 and 500Myr for A- type stars, see Figure 6 of 
ISiegler et al. I [20071 '). In general the levels of excess emission 
are also smaller around lower mass objects (<2 times the 
photospheric flux) apart from around some of the youngest 
sources. These results can be interpreted within the frame- 
work of_the evolutionary models of Kenyon and Bromley 
(2005; 20061 '). These suggest that planetesimals take longer 
to form at 10-30au around A- type stars than at ~lau around 
G-type stars (assuming thermal equilibrium, an excess de- 
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tected at ^24^im implies a temperature of 100-150K, trans- 
lating to dust orbital radii of 3-30au around A and early 
F-type stars and 0.5-3au around solar-type stars). Thus 
there are copious dust producing collisional events for 10- 
lOOMyr during planetesimal accretion around A-type stars. 
Around solar-type stars lOOkm-sized embryos may be com- 
plete within a few Myr, and thus subsequent observed debris 
is likely to have been produced in recent catastrophic col- 
lisions li ke the impac t that formed the Earth-Moon system 
(see e.g. ICanu~pll2004l and references therein). These massive 
collisions are seen in as peaks in the excess emission above a 
lower backgr ound of emission in simula tions of planetesimal 
discs (see e.g. iKenvon fc B romley 2005). Other authors have 
found these collisions could ignite collisional avalanches in 
the disc, which are brighter in discs with high optical depth s 
(large populations of small grains, iGrigorieva et al.l 120071) . 
and could lead to longer-lived increased emission. 

The observed peak in debris disc emission for A 
and early F-type stars (10-20Myr) is an important epoch 
for terrestrial planet formation, as models and radio- 
metric dating suggest that ter restrial planets reach their 
final mass with in 10-30Mvr ilWetherill fc Steward Il993t 
lYin et all |2002| ; IKenvon fc Bromlevl I2006T I. Disc emission 
from this epoch could therefore provide a direct tracer of 
the frequency of terrestrial planet-forming collisions. Cluster 
samples at this age have already been the target of Spitzer 
campaigns. Clusters with published 24/^m Spitzer photom- 
etry include: the 16-17 Myr S co-Cen association (~40 stars 
across F and G spectral types. [Chen et al. 2005); the 27Myr 
old I C 4665 cluster (73 stars, A to mid-K typ es Smith Tet al.l 
l201ll ); NGC 2547 at 35-40 Myrs old (jNavlor fc Jeffriesll2006l ; 
photometry for ~70 sources); NGC 2232 at 25Myr (209 
sour ces B-M type but only 38 have published 24/xm ex- 
cess, [Currie^etj2[2008|); and h and \ Persei at 13 Myr (616 
sources, spectral t ypes B-G but uncertain membership lists, 
ICurrie et al.ll200"8h . It is difficult to obtain sufficient statis- 
tics from individual clusters to examine model predictions 
(either due to uncertain memberships or issues with obtain- 
ing sufficient 24pim photometry). Combining open cluster 
data is the primary method of studying factors such as the 
influence of stellar mass on the incidence of debris discs at 
the epoch of terrestrial planet formation. As open clusters 
provide a homogeneous, chemically uniform coeval popula- 
tion, studies of multiple open clusters allow tests of how en- 
vironmental factors influence debris emission, and the planet 
forming processes that create it. 

NGC 1960 (M36) is a young open cluster , with an age 
deter mined by isochrone fitting of 16l5°Myr l|Sanner et al.l 
2000). The large error bars on this value are a reflection 
of the difficulty of determining the age of this cluster as 
it has no known significantly evolved stars. Other recent 
age estimates in the li terature incl ude an age o f ~20Myr 
|Mavne fc Navlorll2008l ') and 25Mvr (iHasanlfe oOS). As these 
are within the errors given bv lSanner et al.l (|2000). we adopt 
an age of 16l,i Myr in this paper. The cluster lies at a 
distance of 1318 ±120pc and has a reddening of Eb-v = 
0.25 ± 0.02mag (|Sanner et al.l 12000^ . This cluster provides 
an ideal complement to the studies of NGC 2232 and h and 
X Persei as it lies between these clusters in age. Combining 
the studies of these clusters will allow us to place constraints 
on the evolution of debris discs within the peak of debris disc 
emission itself. It is also worth noting that this cluster can 



be viewed as an older analogue of the Orion Nebula Cluster 
(ONC). Integrating a Kroupa mass function l|Kroupall200"lh 
and noting that there are ~30 objects with masses of 3-7M0 
(using a distance modul us of (m — M)p = 1 0.6±0.2mag and 
extinction of Ay = 0.8. ISanner et al1l2000l . and finding ob- 
jects with absolute magnitudes ~-1.4-0.79, see target selec- 
tion in next section) there are ~445 stars with masses 0.5- 
7M@ and ~1680 stars with masses 0.1-7Mq. These num- 
bers are comparable with the ONC. We would also expect 
14 stars more massive than 7M© and find 16 stars of this 
size in NGC 1960 (see target selection in the next section), 
which is also comparable to the ONC. 

In this paper we present Spitzer archive data of NGC 
1960. We use IRAC and MIPS 24^m photometry to identify 
excess emission around cluster members based on Ks — [24] 
vs V — Ks colours and SED fitting. We discuss how levels 
of excess compare with other cluster samples at this epoch, 
and how these observations constrain models of terrestrial 
planet formation. 



2 NGC 1960 TARGETS 

To identify bright target members of t h e clus ter we used 
the sample presented by ISanner et alj (|2000h and adopt 
the same proper motion and colour membership criteria 
presented in that paper. From their initial sample of 864 
stars with B and V band photometry, 121 sources were 
found to have proper motions consistent with cluster mem- 
bership (within 3a of the central cluster proper motion 
2.9 ± 2.7mas/yr in RA, —8.0 ± 2.5mas/yr in Dec). Proper 
motion detections were considered accurate for sources with 
V < 14mag, and thus we adopt this brightness as a cut- 
off for membership from this catalogue (V = 14 trans- 
lates to a spectral type of ~G0 for this cluster). As a 
further test of cluster membership the B — V versus V 
colour-magnitu de of the targets is compared to a 16Myr-old 
isochrone from lSiess et alj (|2000r i with E(B-V) = 25ma g 
(as determined by isochrone fitting in ISanner et al.ll2000l 1. 
The colour-magnitude diagram is shown in Figure [T] This 
led to the rejection of 7 targets with V < 14 and a final 
selection of 63 targets. We extract Ks magnitudes for th e 
targets from the 2MASS catalogue l|Skrutskie et al.ll2006r i. 
This source list is given in Table [TJ 

A set of lower mass members of NGC 1960 was assem- 
bled from spectroscopy obtained at the William Herschel 
4.2-m telescope. A separate paper is being prepared using 
these data, but the details relevant to this paper are given 
below. Observations were performed on 2009 November 25- 
26 and 2010 November 12-14 using the AF2/WYFFOS mul- 
tifibre spectrograph, with an echelle grating that provided 
a spectral coverage of 420A centred at 6600A and a resolv- 
ing power of 11000. Targets with 14 < V < 18.5 were se- 
lected from a V versus V — I colour-magnitude diagram, 
constructed using a BVI photometric survey performed at 
the Isaac Newton Telescope (N. Mayne private communica- 
tion). We took spectra of 472 unique objects (64 were ob- 
served twice or more) along with a suite of standard stars. 

The data were reduced and analysed as described in 
Ijeffries et al.l (|2009f ). Bespoke software was used to Hatfield, 
extract and wavelength calibrate the spectra before cross- 
correlation with radial velocity (RV) standards using the 
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wavelength range 6370-6550A. We obtained good radial ve- 
locities for 358 unique targets - the precision varied with 
signal-to-noise ratio and rotational broadening. In a clus- 
ter as young as NGC 1960, the presence of measurable Li 
in the atmospheres of its cool stars is expected. Indeed, Li 
is depleted so rapidly in low-mass stars that its presence 
in K or M-type stars (V — Ks > 2.0) is a v ery strong in- 
dication of youth (< 100 Mvr, I Jeffrie^ 120061 1 ■ The spectra 
were examined for the presence of a lithium resonance line 
at 6708 A. After correcting all the results to a common helio- 
centric frame, it was clear that the Li-rich objects defined a 
kinematic cluster within the data. After an iterative clipping 
algorithm was applied we found that the Li-rich objects had 
a mean heliocentric RV of — 4.50±0.14kms _1 , with an addi- 
tional estimated zero point uncertainty of about 0.5 kms . 
An analysis of stars with two or more spectra showed that 
our RV precision estimates are realistic. 

To select members we chose to include all those Li-rich 
stars which have a heliocentric RV within 2 error bars of 
the cluster mean, where this error included an additional 
systematic lkms - repre senting a typical in ternal RV dis- 
persion for a cluster (see Ijeffries et all 120091 ). This gave a 
list of 69 candidate members with a median RV precision of 
1.3 kms -1 (targets are listed in Tabled- To estimate the 
level of any possible contamination, we make a similar se- 
lection from the Li-rich stars, but using mean RVs that are 
displaced by ±15kms _1 from the actual cluster mean. This 
resulted in the selection of ~ 6 objects which, despite their 
discrepant RVs, could still be close-binary cluster members. 
We conclude that this sample of low-mass cluster candidates 
has < 10 per cent contamination. 

For both the bright and lower mass samples tempera- 
tures are estimated using the V — Ks colours of the targets 
using the relation 



5040/T cff = 0.555 + 0.195(V - K s ) + 0.013(V 

ijAlonso et al.l 1 19961). Colours were fir st dereddened using 
E(B - V) = 0.2 5 dSanner et al.l l2OO0h and conversions to 
E{V - K s ) from lRieke fc Lebofskvl l| 19851 ). For some of the 
higher mass targets (spectral type earlier than F0) a relation 
in B — V is more often favoured to approximate tempera- 
ture. Using the relation B-V = -3.684 log T cff + 14.551 we 
find temperatures that can be lower by up to 10% for the 
highest mass targets. In this paper temperature is only used 
to determine colour corrections on the photometry. Within 
the ~ 6000 — 11000K temperature range of the higher mass 
targets, this translates to a maximum difference in aperture 
correction factor of <0.5%. 



3 SPITZER DATA 

Data were sourced from the Spitzer Heritage Archive (S HA). 
Data were taken with t he IRAC |Fazio et al.1 l2004h and 
MIPS |Rieke et al.l 12004 ) instruments under Spitzer Pro- 
gramme P50265. IRAC data were obtained on 23rd October 
2008 in high dynamic range mode using a 6x8 mapping ar- 
ray. The AOR consists of 260"map steps, array orientation 
and a 5-point Gaussian dither with small scale factor. A 
frame time of 12 second exposures gave a 5a point-source 
sensitivity of 17.4, 16.6, 15.7 and 14.8 magnitudes in the four 
IRAC bands (3.6/im, 4.5/xm, 5.8/xm and 8/xm) respectively. 



Table 1. Target sources in NGC 1960 from lSanner et al.l te00(t) . 
Magnitudes and colours corrected for redd ening (using A v = 
0.775, E(B - V) = 0.25; ISanner et alj|200d and E(V - K s ) = 
0.686 using the relations from lRieke*fc Lebofskvl Il985l ). Only a 
few rows are shown here to illustrate the content, the full table is 
available online. 



ID 


RA 


Dec 


V 


(B-V) 


(V-JsTs)q 


T e g, K 


1 


84.06585 


34.1436 


8.052 


-0.216 


-0.681 


11769 


2 


84.09608 


34.1758 


8.133 


-0.232 


-0.696 


11842 


3 


84.16359 


34.0639 


8.178 


-0.199 


-0.666 


11696 


4 


84.17627 


34.2017 


8.275 


-0.144 


-0.137 


9535 


5 


84.13335 


34.1798 


8.297 


-0.222 


-0.698 


11852 



Table 2. The lower mass target sources in NGC 1960. Magnitudes 
and colour s corrected for red dening (using A v = 0.775, E(B — 
V) = 0.25; ISanner et aI.ll20od and E( V - I) = 0.400 and E(V - 
Ks) = 0.686 using the relations from iRie kc & Lcbofsky 19851 V 
Only a few rows are shown here to illustrate the content, the full 
table is available online. 



ID 


RA 


Dec 


V 


(V-I)o 


(V-K s ) 


T cff , K 


1.136 


84.06671 


34.0180 


14.378 


0.675 


1.725 


5419 


1.138 


84.04419 


33.8800 


14.050 


0.426 


1.033 


6542 


1.151 


83.91997 


33.9586 


14.155 


0.449 


1.273 


6114 


1.223 


84.09403 


33.9216 


14.999 


0.664 


1.788 


5332 


1.371 


84.12115 


34.0361 


16.086 


0.963 


2.438 


4550 



MIPS data were obtained on 1st November 2008 in 
scan mapping mode centered at 5h36ml2s +34d08m24s. A 
medium scan rate, scan leg length of 0.5° and 11 scan legs 
were used to cover an area of 38'x55'. Five repetitions of 
the AOR were obtained, giving a 5a point-source sensitivity 
of 11.4 magnitudes at 24/im. 

The data were extracted as BCD (basic calibrated 
data) files from the archive. These data are individually 
flux-calibrated array im ages. The Spitzer Science Cen- 
ter MOPEX package |Makovoz fc Marleaul [2005) was 
used to produce the final mosaics. We used standard 
MOPEX modules. The individual 24iim MIPS frames 
were flat-fielded using the flatfield module in MOPEX. 
Overlap correction was determined using the default 
settings in the overlap module and the final image mo- 
saic consisting of all five repetitions of the AOR was 
constructed using the mosaic module. Mosaics were cre- 
ated for each of the IRAC channels using the overlap 
and mosaic modules in MOPEX w ith default settings. 
For d etails of these modules see M akovoz fc Marleaul 
(|2005l) or the on-line M OPEX user's guide at 
http: / /ssc. spitzer. caltech.edu/dataanalysistools/tools/ 
mopex/mopexusersguide/. Photometry was extracted using 
the APEX package from MOPEX. For the IRAC channels 
the PSF is undersampled and thus photometry was ex- 
tracted in circular apertures of radius 3 pixels (~3'/6) with 
background determined in an annulus of inner radius 12 
pixels and outer radius 20 pixels (~14'.'4-24"). Apertures 
were centered on the location of each source as listed in 
Tables Q] and H Of the 132 targets in the full sample, 4 (2, 
4, 2) fell outside the channel 1 (2, 3, 4) imaging mosaics. 
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Figure 1. A V vs B — V colour- magnitude diagram of bright 
candidate members of NGC 1960 from the sample observed by 
ISanner et al,l | |200Gi ) . Only targets with proper motions consistent 
with cluster membership are shown. Proper motion is accurate 
for targets with V < 14, which we adopt as a brightness cut-off 
in our target sample to a void contam i nation by field sources. A 
16Myr-old isochrone from Siess et al. is ovcrplottcd with 

conversion following lKenvon"*fc*^artmannl 1 19951 ) adjusted for a 
distance of 1300pc, E(B — V)=0.25. Targets marked by plus signs 
have colours inconsistent with cluster membership. Crosses mark 
the sources with colours consistent with cluster membership. 



Array location-dependent corrections were applied to the 
IRAC photometry. Correction images available on-line at 
http: / / ssc.spitzer.caltech.edu /irac / calibrationfiles /locationcolor / 
were mosaiced in the same way as the data frames to pro- 
duce a correction mosaic. The aperture photometry of a 
source was then multiplied by the value of the correction 
mosaic at the center of the source image. Aperture cor- 
rections were taken from tabulated values listed in the 
IRAC data handbook. Colour corrections were applied 
by interpolation from tabulated values using the source 
effective temperatures as listed in Tables [T] and [2] We used 
the tabulated values in the IRAC data handbook to convert 
the flux of the targets in Jy to magnitudes. Specifically, the 
zero-points used were 280. 9Jy at 3.6/im, 179. 7Jy at 4.5pm, 
115. OJy at 5.8/im and 64.1Jy a t 8.0/im. Absolute calibration 
of IRAC is stable to 1-3% l|Reach et all 120051 ). We add 
a 3% calibration uncertainty in quadrature to statistical 
background errors from pixel to pixel variation determined 
in the aperture module to give a final error on the IRAC 
photometry. The phot ometry was d e -redd ened using the 
values determined by ISanner et all J2000l) converted to 
IRAC wavelengths using relations of iFlahertv et al.l l|2007t ). 
The final photometry is listed in Tables [3] and [4] 

In the MIPS data the PSF is not undersampled. We 
used the APEX PRF(Point Response Function) fitting mod- 
ule to determine a PRF model for the final mosiac using the 
brightest and cleanest stars (no near neighbours, no bad pix- 
els) in our target list to determine a PRF model using the 
prLestimate module. This model was used to fit the target 
stars to determine source position (from input lists taken 
from Tables [T] and and flux in the PRF photometry mod- 
ule in APEX. If the source could not be well fit (according to 
a x 2 analysis) with either a single or multiple point sources 
(active deblend), then the PRF photometry was determined 
to have failed. For sources which were not fit in the PRF 
photometry, or where the PRF fits were sufficiently distant 



from the input source position that the detected source was 
unlikely to be our target (>l'.'225 corresponding to 1/5 of 
the FWHM of the average PSF for all sources) then the aper- 
ture photometry was adopted instead. Apertures of radius 
2.6 pixels (~6'.'37) were used with annuli of inner radius 8.16 
pixels and outer radius 13.06 pixels (20"- 32") to determine 
the background. 

Some of our targets are close enough to other objects 
for contamination of the aperture photometry (not PRF 
photometry) to be an issue. To determine a correction for 
contaminating sources in the aperture photometry we used 
the APEX module in MOPEX to create a list of all source 
detections (> 3a) in the MIPS 24/im mosaic. We identi- 
fied the sources that were most isolated from other detec- 
tions (no other detection within 10 pixels, 25" of the source) 
and grouped them by brightness (PRF determined flux). 
We then used the aperture module to place apertures (of 
the same size, 6'.' 37 radius, used for the target sources) at 
the source and then at increasing distance from the source 
(in a direction away from other detections) to determine 
the flux that would fall into an aperture at a set distance 
from the source. The relationship between the flux and dis- 
tance from the source was found to be well fit by a semi- 
Gaussian profile with FWHM 10'.'5. Flux from the target 
was <5% of the flux measured in an on-source aperture at 
>10 / /9 and <1% of the on-source flux at >13'.'5. For each 
of our target sources the area within 20"was checked for 
detected sources in the 24pm image, and any possible con- 
tamination calculated by multiplying the aperture flux by 
the value of the Gaussian function described above at the 
appropriate distance. In order to account for sources that 
might be below the 3er detection limit in the 24pm image, 
but that could be contaminating the target photometry, we 
also considered detections recorded by the APEX module 
in the IRAC 8pm mosaic. The relation between 8pm and 
24pm flux was calculated by comparing the values deter- 
mined for sources detected in both images. The median and 
median absolute deviation for the ratio of measured fluxes 
was i 7 8 Mm /i ? 24(jm = 4.57 ± 3.56. Any sources within 20"of 
our target detected at 24pm had their contamination calcu- 
lated, the contaminating source was then removed from the 
list of detected sources at 8pm and then the residual 8pm 
detection list was checked for sources within 20"of the tar- 
get. The flux of these sources was divided by 4.57 to get the 
24pm flux and the contamination in the target aperture cal- 
culated as above. The median and mean absolute deviation 
of the contamination of the targets is 7% ± 6% of the aper- 
ture flux. Sources with high levels of contamination (>13%, 
>lcr higher than average) or with other issues that could af- 
fect their photometry (e.g. high number of bad pixels, close 
to the edge of the mosaic) are noted in Table [3] PRF pho- 
tometry should mitigate against the effects of source con- 
tamination through the active deblend algorithm, and so we 
do not include a contamination correction for sources with 
PRF photometry. We used the above steps to determine the 
levels of contamination that would have been seen in aper- 
ture photometry for those sources with PRF photometry. 
This was used as a method of checking that there could be 
no sources in the IRAC 8pm mosaic that might have been 
blended with the target in the lower resolution 24/im mo- 
saic. In all cases the possible contamination was ^ 1% of the 
PRF photometry. 
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Table 3. Spitzcr IRAC and MIPS 24/xm photometry for on bright members of NGC 1960. Only a few rows are shown here to illustrate 
the content, the full table is available online. 



Source 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 


K s - [24] 


-F24/.Fphot 


Comments 


1 


8.909 [0.033] 


8.844 [0.033] 


8.876 [0.033] 


8.966 [0.033] 


8.918 [0.046] 


-0.183 


0.906 




2 


8.985 [0.033] 


8.953 [0.033] 


8.990 [0.033] 


9.061 [0.033] 


8.784 [0.045] 


0.048 


1.122 




3 


8.968 [0.033] 


8.934 [0.033] 


8.993 [0.033] 


9.034 [0.033] 


9.129 [0.048] 


-0.283 


0.826 




4 


8.106 [0.033] 


7.709 [0.033] 


7.549 [0.033] 


7.310 [0.033] 


6.262 [0.043] 


2.152 


7.424 


XS 


5 


9.154 [0.033] 


9.110 [0.033] 


9.147 [0.033] 


9.221 [0.033] 


9.259 [0.048] 


-0.261 


0.845 




33 


11.140 [0.033] 


11.122 [0.033] 


11.186 [0.034] 


11.254 [0.038] 


>11.453 


<-0.410 


<0.423 


Contam. 



Notes: Stars which had contaminating flux of more than 13% subtracted from the aperture photometry are marked as "Contam" in the 
comments column. None of these targets h ave been identified as a having confirm 24/im excess. Sources with Kg — [24] more than 3cr 
higher than the photospheric relation from IPIavchan et al ] l|2009h are marked as excess sources by the note "XS" in the comments 
column. The error on the photospheric relation includes differences from different photospheric predictions, variation in the colours 
observed for sources without excess and statistical error on the 24/xm photometry added in quadrature as described in the text. 



Table 4. Spitzer IRAC and MIPS 24/^m photometry for low-mass members of NGC 1960. Only a few rows are shown here to illustrate 
the content, the full table is available online. 



Source 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 


K s - [24] 




Comments 


1.136 


12.621 [0.033] 


12.645 [0.033] 


12.679 [0.039] 


12.726 [0.055] 


>1 1.940 


<0.715 


<1.890 




1.138 




12.935 [0.034] 




13.312 [0.119] 


11.886 [0.300] 


1.134 


2.851 


XS 


1.151 


12.809 [0.033] 


12.810 [0.034] 


12.832 [0.040] 


12.838 [0.084] 


11.385 [0.166] 


1.499 


3.958 


XS 


1.223 


13.254 [0.034] 


13.272 [0.034] 


13.284 [0.050] 


13.508 [0.156] 


>12.007 


<1.206 


<2.964 




1.371 


13.580 [0.034] 


13.486 [0.035] 


13.626 [0.058] 


13.828 [0.214] 


>1 1.930 


<1.720 


<4.623 





Notes: See notes on Table [3] for description of comments column. 



After subtracting the contaminating flux from the MIPS 
24/im photometry an aperture correction of 0.433 magni- 
tudes determined from bright targets was applied. A colour- 
correction was applied for each target by interpolation from 
values tabulated in the MIPS data handbook. We used the 
zero-point of 7.14Jy listed in the handbook to convert the 
fluxes into magnitu des. Photometric errors for MIPS 24/im 
observations is 4% |Engelbracht et alj|2007f ). We added this 
in quadrature to statistical errors returned from the aper- 
ture/PRF photometry APEX modules arising from pixel to 
pixel variations to give a final error on the flux. The final 
MIPS photometry is listed in Tables and H Of the 132 
targets in the sample, 2 fall outside the 24/im mosaic (2.195 
and 2.277). These targets also fall outside the field of view 
in the IRAC images, and therefore they are not included in 
the table. 



4 DETECTION OF EXCESS EMISSION 
4.1 Near-infrared excess 

The sensitivity of the IRAC data allows the detection of pho- 
tospheric emission down to the levels of photospheric emis- 
sion of M dwarfs at the distance of NGC 1960 (~1300pc). 
Thus all the targets which fall within the IRAC fields of 
view are detected with good signal-to-noise (>5). 

To determine if there is any excess emission in the IRAC 
photometry from the target stars we plot the Ks — [3.6], 



K s - [4.5], K s - [5.8] and K s - [8.0] colours against 
V — Ks- This method has been used in other clusters to 
search for near-infrared excess (see e .g. iGorlova et al. I l2007l ; 
ICurrie et al]|2008l : ISmith et al.ll201lf ). These plots are shown 
in Figure [2] 

The plots show that the majority of targets have colours 
consistent with photospheric emission. Crosses mark the tar- 
gets in Table [T] and diamonds mark the sources listed in 
Table [5] Errors on the Ks— IRAC colours are taken from 
the pixel-pixel variation in the background apertures in the 
IRAC photometry plus the 3% maximum calibra tion er- 
ror typical for IRAC photometry (|Reach et al.ll2005h , added 
in quadrature to the errors on the Ks photometry from 
2MASS. In order to identify sources with unusual colours, 
we adopt the following procedure. For each target, we iden- 
tify other targets in the sample with similar V — Ks colours 
(within S(V - K s ) < 0.25). The mean 7v" s -IRAC colour for 
these targets is calculated, together with the error (standard 
deviation of the Ks— IRAC colours added in quadrature to 
the mean error on the colour calculated as described above). 
Any targets with colours differing from the mean of other 
targets with similar V — Ks by more than 3 times this error 
are identified and removed from the source list. This pro- 
cess is repeated until no sources with colours different from 
the average by more than 3 times the error are found. The 
outlier sources are labelled on the figures. 

From the IRAC colours, we identify bright targets 4 and 
8 as having colours consistent with excess emission in all 4 



6 R. Smith and R. D. Jeffries 




-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

K 5 - [5.8] K s - [8.0] 



Figure 2. The IRAC colour-colour plots of the targets in the cluster NGC 1960. Errors include pixel- pixel variation and a 3% calibration 
error on the IRAC data, added in quadrature to the uncertainty on the 2MASS Ks band photometry. These plots show colours after 
dereddening (see text for details). Sources with colours indicative of a possible excess are labelled and discussed further in the text. 



IRAC bands, bright target 14 as having colours consistent 
with excess emission in bands [4.5], [5.8] and [8.0], and lower 
mass stars 4_336 and 4_536 as having colours consistent with 
excess in the two longest IRAC wavebands. None of these 
targets are highlighted in the 2MASS catalogue as having 
uncertain Ks photometry. 



Bright target 16 has colours that appear to be consis- 
tent with excess emission in the two longest IRAC bands 
(K s - [5.8] = 0.133 ± 0.038 and K s - [8.0] = 0.249 ± 0.039), 
although in the cr-clipping algorithm described above these 
colours are at the 2.3<r and 2.1<r levels of significance. Con- 
structing an SED for this object using a Kurucz model pro- 
file for the appropriate temperature (9900K, see Table[l} and 
scaling to a best fit for the 2MASS JHKs photometry we 
find the IRAC channel 3 and 4 photometry to be consistent 
with photospheric emission within 2 and 3a respectively. 
Fitting a blackbody to the IRAC and MIPS photometry of 
the source shows that any excess must be at a maximum 
temperature of 600K (in order to not exceed limits from 
the IRAC photometry). This corresponds to a minimum ra- 
dial location of 0.67au assuming the emitting grains behave 
like blackbodies. However, as the IRAC photometry is con- 
sistent with a photospheric emission within the errors we 
cannot confirm that this target exhibits true near-infrared 
excess. 



4-1.1 Contamination by other sources 

The target 4_536 lies very close to the bright target 8. In 
the IRAC channel 4 data (8pm), the emission from bright 
target 8 extends to the location of 4_536, and so the target 
may have contaminated flux. This is also true in the IRAC 
channel 3 image, although the emission from bright target 
8 is less extended in this image. As the resolution is higher 
for the lower wavelength channels, this is not the case in the 
channel 1 and channel 2 images. Thus we cannot confirm 
excess for this target in the IRAC bands. This target is not 
detected at the 3<r level in the MIPS 24/^m imaging, and so 
we have no other photometric measurements to further test 
whether this target exhibits true excess emission. 

The target 4_336 lies very close to a bright infrared 
source (detected in the IRAC images at 84.02458 degrees 
in RA, 34.103 degrees in Dec). This source has been identi- 
fied as IRAS 05327+3404, a probable pre-main sequence K2 
class I/II star with a large optically thick circumstellar disc 
producing a la rge infrared excess (n ot thought to be related 
to NGC 1960, iMaenier et "all 1 19991 ). In the high resolution 
data (IRAC channels 1 and 2) this source can be seen as 
separate from the target 4_336. However, in the IRAC chan- 
nels 3 and 4 and the MIPS 24/mi images, the bright source 
IRAS 05327+3404 extends to cover the location of the tar- 
get 4_336. We can therefore only provide upper limits to the 
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flux for this target in IRAC channels 3 and 4, and in the 
MIPS 24/xm photometry. 



4-1-2 True excess 

The bright targets 4, 8 and 14 (BD +34 1113, BD +34 
1110, and BD +34 1098 respectively) are amongst the high- 
est mass stars in our sample, as indicated by their colour 
in V — Ks and temperature inferred from B — V . These 
sources show evidence for excess emission in the near and 
mid- infrared (see Tableland Figure [2]). We checked these 
sources did not have significant contamination using the 
method described in section 2. We construct SEDs for these 
objects using Kurucz model profiles for the appropriate tem- 
perature for each source (see Table [TJ scaled to a best fit to 
the JHKs 2MASS photometry using a x 2 analysis. These 
plots (Figure [3j show that the near- infrared slopes are not 
consistent with photospheric emission alone. We can fit the 
IRAC and MIPS photometry by adding a single temperature 
blackbody to the photospheric emission for bright sources 4 
and 8. These fits are at temperatures of 800K and 900K re- 
spectively (equivalent to dust lying at radii of 0.73au and 
0.36au respectively assuming the material emits as a black- 
body). For bright target 14 the IRAC and MIPS photometry 
cannot be fit by a single temperature blackbody, and so we 
adopt two blackbody contributions at temperatures of 600K 
and 180K to fit the emission (0.73 and 8.06au). It should be 
noted that this illustrative fit is non-unique and that longer 
wavelength photometry would be required to determine con- 
straints on the cooler blackbody temperature in particular. 
In all three plots the blackbody contributions are shown by 
a dotted line and the total (blackbody and photospheric 
emission) is shown by a dashed line (Figure [3]). 

For all three targets with near-infrared excess we must 
ask if this emission is evidence for a remnant primordial 
disc. At an age of 16Myr for the cluster such primordial 
discs would be rare as they are expected to dissipate on 
timescales of a few Myr (see IWvattl 120081 and references 
therein). The fractional luminosity of these discs is much 
lower than we would expect for primordial disc emission. 
Adopting the blackbody fits to the emission shown in Figure 
[3] we find that the excess emission has fractional luminos- 
ity La/L* = 3.9 x 10~ 3 for both bright star 4 and 8, and 
that the two components of the excess around bright star 
14 have fractional luminosities of Ld/Li, — 1.3 x 10 -3 and 
2.7 x 10 -3 (for components at 600K and 180K respectively). 
These levels are consistent with L c \/L* < 10~ 2 which con- 
venti onally defines a debri s disc as opposed to a primordial 
disc (|Lagrange et alJl200Ch . These targets join a small but 
growing pop ulation of A stars with debris disc s within lOau 
(e.g C Lep. iMoerchen et al.ll2007l , HD172555 iRebull et al.1 
2008). Bright star 14 is one of a growing number of systems 
with eviden ce for multiple com ponent discs. 11 stars were 
identified bv lChen et alj (|2009l l as having probable multiple 
component discs, which have been interpreted as possible 
solar-syst em analogues (see also the /3 Pic association mem- 
ber T) Tel, ISmith et alj|2009f j. 
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Figure 4. The spectral slopes for all targets in our study. Sources 
without significant detections at 24/im are shown as upper limits 
to the spectral slope based on the 3<r upper limits to their flux. 
Regions of the spectral slope parameter space which could indi- 
cate primordial (optically thick) or transition discs are labelled. 
Most targets have spectral slopes consistent with no excess or 
optically thin emission consistent with a debris disc. See text for 
discussion. 



4-1.3 Limits on near-infrared excess and primordial 
emission 

As a further confirmation of the evolved nature of the 
targets observed in th is study, we follow the example of 
IMuzerolle et~ai1 l|20ld ) and consider the IRAC/MIPS SED 
slope for each of our objects as a method of identifying pos- 
sible primordial/transition disc candidates. We consider the 
value of a = dlog vF/dlog A over IRAC channels 1 and 3 
(so for A = 3.6/im and 5.8/im) and over IRAC channel 4 
and MIPS (A = S.Ofim and 23.7/xm). We plot these spectral 
slopes for all stars in our study in Figure 2J and plot upper 
limits to the spectral slope as-24 for sources without sig- 
nificant photometry a t 24^im. According to the metrics of 
IMuzerolle etlll (|2010T ). stars with 03.6-5.8 > —1.8 are pos- 
sible primordial disc hosts. Stars with 03.6-5.8 < — 1.8 and 
«8-24 > are possible transition disc candidates, and may 
have optically thick discs at >lau with depleted/optically 
thin inner holes. Target 4_336 appears to lie in the primor- 
dial disc region, but this object's flux is highly contaminated 
by a known primordial disc candidate as discussed in Sec- 
tion 4.1.1, and therefore this is not a true detection of a 
primordial disc. Bright target 4 lies close to the boundary 
for possible primordial emission. Bright target 14 lies in the 
region of possible transition disc emission. However, for both 
these targets we have shown that it is possible to fit their 
emission with optically thin (Lir/L* < 10~ 2 ) blackbody 
emission. Thus we find the emission for both these targets 
is consistent with debris disc emission. The remaining tar- 
gets in our study all have 0:3.6-5.8 < —1.8 and O8-24 < 0, 
consistent with no excess or optically thin emission. 

For those targets with no evidence of significant near- 
infrared excess, we determine the limits we can place on 
excess emission in the IRAC photometric bands. Exclud- 
ing the sources discussed above, we determine an empirical 
photospheric locus and uncertainty in Ks — IRAC as a func- 
tion of V — Ks- For each value of V — Ks, we determine the 
mean and standard deviation of the Ks — IRAC colours 
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Figure 3. The SED of bright targets with significant excess 
emission in IRAC channels. Crosses mark Hipparcos V band 
and 2MASS JHKg photometry. Diamonds mark IRAC and 
MIPS 24/j.m photometry. A Kurucz profile of appropriate tem- 
perature for each source (see Table [TJ is shown with a solid 
line. Dotted lines mark blackbody fits to the excess emission, 
with the total (blackbody plus photospheric) emission indi- 
cated by a dashed line. 
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Figure 5. The Kg — [24] vs V — Kg colour-colour diagram for sources with > 3cr detections in the MIPS mosaic. Grey symbols mark 
uppe r limits for sources with < 3tr detections. Diamonds mark sources from the lower mass sample, crosses mark additional sources 
from ISanner et al.l (l2000h . Errors are the statistical and calibration e rror o n [24] added in quadrature to the error on the 2MASS Kg 
magnitude. The predicted photospheric colours from [Plavchan et al. (2009) are shown with a solid black line. For most of the spectral 
range under consideration, only sources with significant excess are detected. 
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of sources with V — Ks within 0.25 of this value. These 
are our empirical photospheric locus and error. We fit a lin- 
ear model to the error to determine the sigma values for 
each colour. We determine the dispersions in K s — IRAC 
colours to be <tk s -[3.6\ = 0.04 + 0.005(U — Ks), o\k s -[4.5] = 
0.04 + 0.01(V - Ks), (T Ks -[B.e\ = 0.04 + 0.02(V - K s ), and 
tk s -[8.o] = 0-05 + 0.04(V — Ks)- Adopting these uncertain- 
ties, we find that for an F5-type star (V — Ks ~ 1) we can 
place the following 3<r limits on excess emission in the near- 
infrared: Fs.e/Fphot < 1.13; F4. 6 /-F P hot < 1.15; F 5 . 8 /F phot < 
1.18; and F s .o/F phot < 1.28. For objects with V - K s = 4 
(~ Ml-type) we place the following 3er limits on excess in 
the near-infrared: F3. 6 /F phot < 1.18; F 4 .5/F phot < 1.25; 
Fs.s/Fphot < 1-39; and F 8 . /F phot < 1.79. 

4.2 24/im excess 

The principal means of identifying debris discs is the de- 
tection of excess emission (above the levels expected from 
the photosphere) in mid-infrared ranges, and no detection 
of such an excess at near-infrared wavelengths (which may 
indicate a primordial or transition disc). The sensitivity of 
the MIPS observations of NGC 1960 allow a detection down 
to ~12th magnitude at 24pm (3cr point source limit). This 
means we can only detect photospheres to A3V (V — K < 
0.10), which means that for many of our targets, we would 
not expect to detect the photosphere. This is reflected in 
the results presented in Tables [3] and [4] Those sources with 
significant detections generally have high values of Ks — [24] 
and -F24/fphot (detected emission at 24/im / expected emis- 
sion from the photosphere). 

In order to determine which targets exhibit 24/im 
excess, we follow the example of recent authors (e.g. 
iRebull et alJl2008l : IStauffer et ai]l20ld ) in using the targets' 
Ks — [24] colours. This m ethod requires a well- defined model 
for photospheric colours. IStauffer et al.1 (|2010h used Spitzer 
observations to determine an empirical relation for Ks — [24] 
from V - K s , defined as K s ~ [24] = 0.042 - 0.053 x (V — 
Ks) + 0. 023 x (V — Ks ) 2 . Sim ilar r elations have been pre - 
sented bv lGorlova et al.1 l|200d ). and lPlavchan et al.1 lj2009h . 
who included M dwa rfs in the spectral range covered. As the 
IStauffer et all (|2010l ) relation is only valid to V — Ks < 3 
and our NGC 1960 sample incl udes sources of lower mass, we 
use the iPlavchan et all (|2009h relation as our photospheric 
model. 

In Figure [S] we show the Ks — [24] vs V — Ks colours of 
the target source s. Crosses are used t o mark the bright tar- 
gets taken from ISanner et all (|200Cf ). and diamonds mark 
the sources identified as lower mass cluster members. Up- 
per limits for targets with signal to noise of less than 5 are 
marked in grey. Over-plo tted as a solid line is t he expected 
photospheric colour from IPlavchan et al.1 (|2009h . We deter- 
mine the dispersion around this relation by considering the 
colours of the bright targets (V — K s < 0.1) for which the 
photosphere should have been detectable. A histogram of 
the Ks — [24] colours of these sources was constructed, and 
a Gaussian function fitted to the left-hand side of this func- 
tion (the targets which are not likely to have excess emis- 
sion). The best fitting Gaussian FWHM was found to be 
0.14. This dispersion includes the errors in Ks — [24] (typi- 
cally around 0.05 for targets with V — Ks < 0.1, see Table 
[3|, and so we adopt this value as a minimum uncertainty on 



the photospheric values of Ks — [24]. We expect the errors 
on the Ks — [24] colours to increase for fainter targets due 
to decreased signal to noise. Using all significant detections 
(> 5a) we examined the errors on the 24/im photometry as 
a function of V — Ks- We found that the effect of source 
colour on the uncertainty was relatively weak. This is be- 
cause for most of the spectral range considered in this study 
we are detecting sources with excess emission only, and so 
the error levels were nearly constant. We determined a lin- 
ear fit to the uncertainty (using a x 2 test) on Ks — [24] and 
added this to the measurement of the dispersion adopted 
as a minimum uncertainty as discussed above. This gave us 
a final uncertainty on Ks — [24] from our photometry of 
0.14 + 0.032(V — Ks)- We also consider th e differences be- 
tween t he photospher i c colo ur m odels from | Plavchan et al.l 
|2009l ). lGorlova et al.1 l|2006l ) and lStauffer et al.1 (^OloT ). For 
each value of V — Ks we interpolate the predicted Ks — [24] 
for each model, and add the standard deviation of the three 
results in quadrature to the error relation already described. 
These errors are shown (at a 3<r level) by a dotted line in 
Figure \5\ 

From this figure we identify 17 of the targets from Table 
[1] and 18 of those from Table [2] as having significant 24/im 
excess. These sources are identified in Tables and [3] by a 
note in the Comments column. 



5 DISCUSSION 

5.1 Relationship with multiplicity 

As the ma j ority of stars are binary or higher-order multiples 
jAbtll 19831 : iDuquennov fc Mavor|[l99ll 'l , the relationship be- 
tween stellar multiplicity and debris d isc parameters ha s 
been the subject of much recent study. ICieza et al.l (|2009l ) 
presented an analysis of IRAC data from several Spitzer 
legacy surveys which showed that for projected separations 
of <40au, systems were half as likely to retain their primor- 
dial discs than systems with larger separations (suggesting a 
disc lifetime of 0.3-0.5 Myr for close bi naries compared to 3 - 
5 Myr around single stars). Conversely iTrilling et al.1 (|2007T l 
found in their study of field stars (with most stars >600 
Myr old) that the binary stars had a higher incidence of de- 
bris discs than single stars. Debris discs were found to be 
more common around binaries with small (<3au) or large 
(>50au) separations than around intermediate separation 
binaries (3-50au). This survey concentrated on detections at 
70/im, and so on discs that were further from their host stars. 
|Plavchan~e t al. (200 §) f o und n o evi dence of the tren d sug- 
gested by ITrilling et al.1 ([20071 ), and iDuchenel (|2010h found 
no significant dependence of debris disc inc i dence with bi- 
narity or binary separation. IStauffer et all (|2010T ) showed 
tentative evidence that in the 100 Myr Blanco 1 cluster bi- 
nary stars have a lower 24/im excess frequency than single 
stars (adopting height above a single star isochrone as a bi- 
narity proxy). They combine their results with data from 
the Pleiades (~ 100 Myr) and NGC 2547 (~ 35 Myr) and 
find an overall chance of 0.05% that the excess around single 
and binary samples is draw n from the same pa rent popula- 
tion (using a K-S test). In ISmith et alj I^Oll! ) we showed 
that there was no evidence for a dependence of debris disc 
emission on stellar multiplicity in the 27 Myr old open clus- 
ter IC4665. We adopt the same analysis here using height 
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Figu re 6. Dependen ce of excess on a proxy for stellar multiplicity. Left: Height above a single star isochrone of 16Myrs is used following 
IStauffer et"aflj2010h . Diamonds mark the targets listed in Table [2] and crosses the targets in Table [l] Stars identified as having excess 
emission (see Tables [3] and |4j are overplotted with large circles. Magnitudes and colours are corrected for reddening and extinction. Right: 
Height above the isochrone is compared with K$ — [24] colour. There is no evidence of a significantly different distribution of colours for 
suspected multiples (V — V prc d > 0.3). 



above a single star isochrone as a proxy for multiplicity and 
compare this to Ks — [24] colour. Th e isochrone wa s tune d 
from a fit to the Pleiades following IStauffer et all (|201CJ| . 
This is shown in the left-hand panel of Figure [6] Bright 
star 4 (HIP 26354) stands out in this plot as having a very 
high V band magnitude for its colour (V — Ks)- This source 
has been identified as having high variability (5V = 0.281 
iLefevre et alJl2009T ) and therefore the target's photometric 
measurements are dependent on the epoch of observation, a 
possible cause of the anomalously high V band magnitude. 

We expect that most targets will have V — Vp, re d ^ 0.75 
(where Vp re d) is the predicted V band magnitude of a star 
of a given V — Ks colour from the single star isochrone). 
Sources with V — V ple a > 0.75 are possible triples or higher- 
order systems. From our sample of 132 targets (from Ta- 
bles [2] and [1} 26 are found to have V — V pfe d > 0.75, giv- 
ing a detection rate of 19lg% higher-order multiples. This 
is somewhat higher than is seen in other studies (~9-1 2%, 
lAbtlll983l ; lAbt fe Lewi 19761 iDuauennov fc Mavoj|l99lf ) , al- 
though not significantly so. This could also be reduced if 
we take into account errors on the 2MASS Ks photometry, 
which means that at the 3a level a minimum of 16 sources 
(12% of the target list) have V — V pIe d > 0.75, which agrees 
with the expected frequencies of higher-order multiplicity. 
In general we do not see any evidence for a relation between 
V — Vp TC d and Ks — [24] (see right-hand panel of Figure HU. 
We separate the sample into sources with V — V pl ed ^ 0.3 
and those with V — Vprcd > 0.3 and compare the Ks — [24] 
colours of the two populations. Using a K-S test on the 
Ks — [24] excess for sources with significant 24/im detections 
(where the excess is defined as (Ks — [24] ) b s — {Ks — [24] ) pro d 
and (Ks — [24]) prc d is determined from the Plavchan model 
for the photospheric colours) we find only a 70% probabil- 
ity that we can reject the null hypothesis that the samples 
are drawn from the same population. We also compare the 
two samples (V — Ks ^ 0.3 and V — Ks > 0.3) includ- 
ing the upper limits. We use the ASURV statistical package 
l|LaVallev et al.|[l990l) to compare the two p opulations using 
several tests (see Feigelson fc Nelsonlll985l for a description 
of the tests). We find that the probability of rejecting the 



null hypothesis (that the samples are drawn from the same 
population) is 0.47 according to the logrank test, and 0.74 
according to the Peto-Prentice test. The differences in the 
probabilities arise from the different treatments of the upper 
limits in the methods. In either case we see no evidence that 
the multiplicity of the stellar system affects the incidence of 
debris disc emission in this cluster. 



5.2 Placing NGC 1960 into context 

The NGC 1960 cluster is at an interesting age for planet 
formation. At 16Myr, any gas giant planets are expected to 
have formed and the circumstellar disc evolved from a gas- 
rich to gas-poor state wit hin ~ 3Myr for high-m ass stars (B, 
A and early F-type, e.g. Hernand ez et al.l 2005h or 5-7Mv r 
for low mass stars (F5 or later, e.g. Hernandez et alj|2007l ). 
Terrestrial planet formation on t he other hand is expec ted 
to continue for up to lOOMyr (e.g. rWeidcnschillingl l 19771 ). A 
peak in 24pm debris emission at 10-20 Myr around A and 
early F-type stars can be interpreted as evidence of delayed 
stirring of the disc by the formation of Pluto-si zed bodies 
at this epoch (see e.g. iKenvon fc Bromlevl |20o3 ) - We now 
explore how NGC 1960 fits in with these models. 



5.2.1 Excess emission at 24fJ,m 

A 'rise and fall' in 24pm excess emission (emission 
above expected photos pheri c levels) was f irst id entified by 
IHernandez et all l|2007l ) and lCurrie et al.l (|2008l ). The peak 
levels of excess emission (defined as emission above that aris- 
ing from the stellar photosphere) and frequency of excess are 
seen to rise from ~5-10 Myrs, before falling after ~ 20 Myrs 
around higher mass (A and early F-type) stars. We examine 
NGC 1960 in the context of this pattern by comparing the 
rate of excess around the higher mass stars to other clusters 
of ages ^5 Myrs. 

In Figurc[7]we show the frequency of excess emission for 
stars of spectral types B0-F5 observed with MIPS. We in- 
clude data for th e following op en clusters: A Orionis (5Myr, 
IHernandez et all 120091 . l201oh : Upper Scorpius (5 Myr, 
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ICarpenter et all 120091 and IChen et al.l 120051); Orion OBlb 
(5 M yr) and Orion OBla (8.5Myr, both iHernandez et all 
2006); Lower Centaurus Cr ux (16Myr) and Upper Cen- 
taurus L upus (17Myr both Ichen et alj l2005h : NGC 2232 
(25Myr, ICurrie et all l20Qg| . note that only 38/209 stars 
have published M IPS 24/xm observatio ns); IC 4665 (27Myr , 
ISmith et al.ll201ll ) ; NGC 2547 (35Myr , IGorlova et al]|2007f ) ; 
IC 2391 (50Mvr, ISiegler et all 12007m ; Blanc o 1 (lOOMyr 
Stauffer et"ai1 l2010l ); the P leiades (115Myr IGorlova etaD 
20061 ); the Hvades (625Mvr, ICieza et alJlioog ); and a collec- 



tion of field stars (|Brvden et alj 20061 ). For all clusters we 
extract the spectral type for each object listed in the above 
papers. Where no spectral type is listed we assign a spec- 
tral type based on the V — Ks colour of the target. The 
relationship was determined from Kurucz model profiles of 
stars at typical temperatures for each spectral type (so F 5 
stars have V — Ks ~ 1, temperatures taken from ICoxl f2000 ) . 
The frequencies are listed in Tabled Clusters with less than 
10 members in the spectral range B0-F5 were excluded from 
Figure [7] but are listed in Tableland included in our exami- 
nation of levels of excess emission below. Stars are identified 
as having excess emission if th ey satisfy the following criteria 
(following ICurrie et al.ll2008f ): 



(K s - [24]) obs - (Ks - [24]) p hot ^ 3a [24] 
and 

(Ks - [24]) obs - (K s - [24]) phot > 0.15, 



(1) 



(2) 



where this last condition takes into account systematic er- 
rors, e.g. from dispersion in Ks — [24] colours. It should 
be noted that all the excess sources identified in NGC 
1960 meet these criteria. (Ks — [24] Y B hot for each target i s 
given by the relation determined bv lPlavchan et all l|2009h . 
as shown in Figure [5] For this spectral range we expect 
Ks— [24] ~ (dereddened). Extinction (Ay) was taken from 
the papers mentioned above, and converted to Ak q and 
according to the rela tions in iRieke fc Lebofskvl (| 19851 ) and 
iFlahertv et al.l (|2007l ) respectively. 

In Figure[7|and Table[S]we give a la confidence limit on 
the frequency of excess. We assign a confidence limit based 
on the t abulated confid ence limits of Poisson statistics pre- 
sented in lGehrelsl <|l986l ). However, for several of the samples 
listed in Table[5]there are sources with 24/im upper limits or 
non-detections, which could potentially be disc hosts. For all 
sources with upper limits to their 24/im photometry, we take 
the 3<r upper limit to [24] and determine if the source could 
have an excess. This is determined by checking if the upper 
limit would be listed as an excess according to equations [1] 
and[2] If the upper limit is below these cut-offs, the source is 
an upper limit with no disc. Otherwise the target is a poten- 
tial host of excess emission. We use the number of sources 
with possible excess emission to define another upper limit 
to the confidence limits on the frequency of excess for each 
cluster. The largest of this upp er limit a nd th at determined 
from the Poisson statistics of iGehrelsl (l986) is presented 
in Table [S] and shown in Figure [7] In this plot we can see 
that the cluster NGC 1960 has a fairly typical frequency 
of excess emission given its age (30l^ 4 %). We can possibly 
see the peak of excess frequency at ~10-20 Myr, however 
the statistics for the younger clusters are quite uncertain. 
As only one young cluster (Upper Sco, ~5 Myr) has a low 
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Figure 7. The frequency of 24/im excess emission of stars of 
spectral types B0-F5 as a function of age. Clusters are listed in 
the text and in Table |51 The datum for NGC 1960 is shown by a 
diamond and is offset in age to 18 Myr to avoid confusion with 
data from the Sco-Cen association. 



frequency of excess, we cannot confirm with this data that 
the "rise and fall" in frequency is statistically significant. 

In addition to considering the overall frequencies of ex- 
cess, we also consider the level of excess emission for each 
target, given as [24] prc d — [24] bs, where [24] pre d is determined 
from predicted (Ks — [24]) p hot for the target (see above para- 
graph) . We show the excesses for stars in the spectral range 
B0-F5 in Figure [8] Here we incl ude observed exces s emis- 
sion from stars in h and \ Persei (jCurrie et al.|[20o3 ) which 
were excluded from the statistical analysis because of un- 
certain membership lists, and j3 Pictoris and the field stars 
which were excluded from Figure [7] because of the low num- 
ber of objects (see Table |S}. Here we again see that NGC 
1960 (shown by diamonds) is comparable to other clusters 
of a similar age. Bright star 14, which has the highest lev- 
els of excess emission in this cluster (and also has excess 
in the near infrared) has one of the highest levels of excess 
emission observed, although its SED shows that the emis- 
sion is consistent with debris disc emission rather than that 
from a primordial disc (see section 4). This figure displays 
more clearly that there is a 'rise and fall' in the upper en- 
velope of 24/im excess emission, as the peak levels of 24/im 
excess reach a maximum at ~10-20Myr before falling off. 
A fall-off proportional to time is e xpected from models of 
the evolution of d ebris discs (e.g. iDominik fc Decm]|2003l ; 
IWvatt et af1l2007l ). 

We show the same plots for lower mass stars of spec- 
tral types F6-K9. In Figure [9] we can see that the frequency 
of 24/im excess falls off roughly in proportion to time af- 
ter ~50Myr. NGC 1960 has a fairly typical excess frequency 
for its age (23 +69 — 9%), although given the large upper 
limit, we have really defined a lower limit to the frequency 
of >23% (excluding the Poisson statistical errors). There 
is some evidence of the 'rise and fall' of debris discs seen 
in higher mass stars, although again the statistics for the 
younger clusters are highly uncertain. In this plot it ap- 
pears that the peak frequencies of debris discs around F6- 
K9 type stars are seen in clusters of ages ~ 30 — 50 Myrs 
with an abrupt fall-off at ~100 Myrs (see also Table[5]). This 
is slower than the expected time taken to form Pluto-mass 
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Table 5. The excess frequencies of open clusters with published 24/jm MIPS photometry. 



Name 


Age 


B0-F5 


F6-K9 


Reference 




Myr 


Excess /Total 


% with excess 


Excess/Total 


% with excess 




A Orionis 


5±1 


16/44 


36 +4S 
dD -9 


22/401 


tr+90 

°-i 


Hernandez et al. (2009, 2010) 


Orion OBlb 


5±1 


11/28 


oq+16 
o»_ 12 


3/5 x 




Hernandez et al. (2006) 


Upper Sco 


5±1 


10/70 




6/23 




Carpenter et al. (2009); Chen et al. (2005) 


Orion OBla 


8.5+.1.5 


12/25 




1/1* 




Hernandez et al. (2006) 


Pic 


12±6 


5/9 x 




2/12 


-,-+29 


Rebull et al. (2008) 


Sco-Cen 


16-17 


13/25 


52+ 19 


3/11 


97+27 
Z '-15 


Chen et al. (2005) 


NGC 1960 




21/71 


30+ 24 


14/61 


23+f 


This paper 


NGC2232 


25±4 


9/23 


oq+26 


4/15 


26+ 21 


Currie et al. (2008) 


IC4665 


27±5 


9/23 


OQ + 18 

° 3 -13 


12/29 


41 +16 
^ 1 -12 


Smith et al. (2011) 


NGC2547 


35±5 


15/33 


45+15 


21/41 


rci+14 
0i -ll 


Gorlova et al. (2007) 


IC2391 


50±5 


2/18 


11 + 15 
iJ --7 


6/13 


46+ 28 
* D -18 


Siegler et al. (2007) 


Blanco 1 


100±20 


4/12 


qq + 26 
°°-16 


1/25 


*±S 


Stauffer et al. (2010) 


Pleiades 


115±20 


7/44 




3/38 




Gorlova et al. (2006) 


Hyades 


625±50 


1/11 


si! 1 


1/67 


it? 


Cieza et al. (2009) 


Field 


4000 a 


0/4* 




6/65 


q+6 
y -4 


Brvden et al. (2006) 



Notes: x Frequency of excess emission in this spectral range highly uncertain due to low number of stars, and so this point is omitted 
from Figure 171 or p9l 

a Average age for all stars in this paper. 
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Figure 8. The 24/xm excess emission of stars of spectral types B0-F5 as a function of age. Clusters are listed in the text. Data for NGC 
1960 are shown by diamonds and offset in age to 18 Myr to avoid confusion with the Sco-Cen association. Grey symbols mark 3cj upper 
limits for sources with low signal to noise 24/^m data. 
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Figure 9. The frequency of 24/^m excess emission of stars of 
spectral types F6-K9 as a function of age. Clusters are listed in 
the text and in Table |51 The datum for NGC 1960 is shown by a 
diamond and offset in age to 18 Myr to avoid confusion with the 
Sco-Cen association. 



pr oto-planets at ~1 au , and in fact according to the models 
of iKenvon fc Bromley! (|2005f) . a later peak in the 24^tm ex- 
cess emission is to be expected for discs at greater distances 
from the star (see Figures 7 and 8 of IKenvon fc Bromlevl 
2005). This suggests that rather than witnessing excesses 
arising from terrestrial planet regions around F6-K9 type 
stars, we may be seeing the inner edge of more distant plan- 
etesimal belts (due to the apparently smooth increase in 
frequencies towards a peak at 30-50 Myrs). Constraints on 
longer wavelength emission for the cl uster samples would 
be required to explore this possibility. Sie gler et al.l (|2007m 
found evidence that for clusters older than lOMyrs the fre- 
quency of 24^tm excess emission fell-off steadily. We do not 
find such a strong relation in our data (we have considered 
a larger set of cluster samples), although it is clear that by 
lOOMyr excess frequencies are very low. These differences are 
likely to arise from different criteria for identifying the pres- 
ence of excess. This plot demonstrates that the frequency of 
debris disc emission is likely to be a function of many fac- 
tors aside from the age of the cluster, however with the small 
numbers of clusters under consideration, an examination of 
these factors is beyond the scope of this paper. 

In Figure [10] we see that the levels of excess emission 
observed around stars in the spectral range F6-K9 in NGC 
1960 are typical for stars of similar ages. We do not see 
evidence for a delayed peak (levels of excess emission be- 
ing higher in clusters ~ 10 — 20 Myr old) in this spectral 
range. This may in part be due to residual primordial disc 
populations at 5Myr around lower mass stars, for exam- 
ple those seen in Upper Sco (sources in this cluster with 
[24l Dre d — [24l bs > 4 were determined to be primordial in 
ICarpenter et all 120091 ). However removal of the primordial 
disc targets would mean that the frequency of true debris 
discs in these clusters is lower, and therefore the evidence 
for a delayed peak would be stronger. In common with the 
drop off in frequency of excess in this spectral range, high 
levels of excess emission are not seen at ages ^ lOOMyrs in 
these clusters. 

Interestingly, for NGC 1960 we see no evidence that 
either the frequency or the levels of excess emission at 



24pm are dependent on spectral type. This is somewhat 
contrary to the expectation that 24pm excess should be 
longer lived around brighter stars, which arises from the dif- 
ference in time taken to form planetesimals at larger radial 
distances around A-type stars (excess detected at ^24/im 
imply a temperature of 100-150K, translating to an off- 
set of 3-30au around A and early F-type stars and 0.5- 
3au around solar- type stars assuming thermal equilibrium). 
Within the context of current planet formation theories 
l|Kenvon fc Bromlevll2005l.l2006r ). we might expect the 24/xm 
emission around lower mass stars of ~16Myr to be the result 
of massive collisions between proto-planetary bodies, and 
we might expect it to occur less frequently than the plan- 
etesimal collisions thought to be responsible for the 24pm 
emission around higher mass stars. In fact, for the majority 
of the clusters listed in Table [5j there is no significant dif- 
ference between the frequency of excess emission observed 
around B0-F5-type stars and that observed around F6-K9- 
type stars. However, for the majority of clusters older than 
>5Myr, the typical levels of excess emission observed around 
higher mass stars are higher than that observed around lower 
mass stars (see Figures 151 and [TO)) . which does agree with the 
model predictions. This does not appear to be the case for 
NGC 1960. The highest observed excess is indeed observed 
around a higher mass star (bright target 14), but in general 
the levels of excess emission observed are similar regard- 
less of V — Ks colour (proxy for spectral type). Our results 
are biased by the fact that only relatively high levels of ex- 
cess will be detected around lower mass stars, and so there 
is probably a population of stars in the range F6-K9 with 
lower levels of excess than shown in Figure [TOl However, as 
we do not see a population of higher mass stars with larger 
excesses (excluding bright star 14), these results might sug- 
gest that there are a large number of lower mass sources in 
the NGC 1960 cluster that are currently undergoing mas- 
sive collisions. The reason why many sources in the cluster 
should experience massive collisions at a similar time is un- 
clear, and if this is the reason for the 24pm excess observed 
around lower mass stars, this presents a challenge to current 
planet formation theories. In Table [6] we show the frequen- 
cies of stars with large 24pm excess ([24] prc d — [24] Q b s > 2) 
for clusters with populations >10 stars in both spectral bins. 
From these statistics we cannot see a significant difference 
between the frequencies of large excess emission observed 
in NGC 1960 and other clusters. We do find some evidence 
that the frequency of stars with large excess emission drops 
at ^25Myr amongst both high and low mass samples, which 
is in agreement with steady-state evolution models. We can- 
not determine if the relationship between the frequency of 
large excesses in the lower mass (F6-K9 type) and higher 
mass (B0-F5 type) spectral bins is significantly different for 
NGC 1960, partly due to the high errors involved (errors on 
the frequencies were determined in the same way as errors 
on the overall excess frequencies in Table [SJ. With the cur- 
rent data, we have no statistically significant evidence that 
NGC 1960 is different from other clusters of similar ages. 

In summary the levels and frequency of excess emission 
seen in NGC 1960 are not significantly different to other 
open clusters of similar ages. The bright target 14 (BD +34 
1098) has one of the highest levels of excess emission ob- 
served in early-type stars, but SED fitting has shown that 
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Figure 10. The excess emission of stars of spectral types F6-K9 as a function of age. Clusters are listed in the text. Data for 

NGC 1960 are shown by diamonds and offset in age to 18Myr to avoid confusion with the Sco-Cen association. Grey symbols mark 3<r 
upper limits for sources with low signal to noise 24/im data. 



Table 6. The frequencies of large excess emission ([24] prct j — [24] t, s > 2) among open clusters with published 24/im MIPS photometry. 



Name Age B0-F5 F6-K9 Reference 

Myr Excess/Total % with excess Excess/Total % with excess 



A Orionis 


5±1 


0/44 


0+ 4 
u -o 


4/401 


111 


Hernandez et al. (2009. 2010) 


Upper Sco 


5±1 


2/70 


3+ 4 


6/23 




Carpenter et al. (2009); Chen et al. (2005) 


Sco-Cen 


16-17 


5/25 


20+^ 


1/11 


9+f 


Chen et al. (2005) 


NGC 1960 


16±l° 


3/71 


4+4 


8/61 


13±t 


This paper 


NGC2232 


25±4 


1/23 


4 +l0 
-3 


0/15 


0+12 
U -0 


Currie et al. (2008) 


IC4665 


27±5 


0/23 


+8 
U -0 


0/29 


0+6 


Smith et al. (2011) 


NGC2547 


35±5 


0/33 


+6 
U -0 


2/41 


°-3 


Gorlova et al. (2007) 


IC2391 


50±5 


0/18 


+10 


0/13 


0+14 
U -0 


Sieeler et al. (2007) 


Blanco 1 


100±20 


0/12 


0+15 


0/25 


0+ 7 


Stauffer et al. (2010) 


Pleiades 


115±20 


0/44 


0+ 4 

U -0 


0/38 


0+5 
U -0 


Gorlova et al. (2006) 


Hyades 


625±50 


0/11 


0+17 


0/67 


0+3 


Cieza et al. (2009) 



Notes: Errors on excess frequencies are determined in the same manner as for Table[5]as described in the text. Only clusters with 24/im 
observations of > 10 stars in each spectral bin have been included in this table. 



this is consistent with debris as opposed to primordial emis- 5.2.2 Excess emission at 8/j,m 
sion. 

A subset of the open clusters discussed above also have pub- 
lished 8pm IRAC photometry. We consider how the fre- 
quency of 8/im excess evolves with time. Frequencies are 
listed by age and spectral bin in Table [7] Errors on the fre- 
quency are derived from upper limits and the Poisson statis- 
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Table 7. The excess frequencies of open clusters with published 8/im IRAC photometry. 



Name 



Age B0-F5 

Myr Excess/Total % with excess 



F6-K9 Reference 
Excess/Total % with excess 



A Orionis 


5±1 


2/44 


5+ 6 
°-3 


15/401 


411 


Orion OBlb 


5±1 


1/29 


q + 34 

J -3 


0/5 x 




Upper Sco 


5±1 


5/68 


7+5 
' -3 


4/45 


9t S 5 


Orion OBla 


8.5±1.5 


2/25 




0/l x 




NGC 1960 


16^ 


3/71 


*±t 


0/61 


<7 


NGC 2232 


25±4 


1/23 


411° 


1/15 


7+15 
'-6 


IC4665 


27±5 


1/23 




0/29 


<10 


NGC 2547 


35±5 


2/33 


6±g 


2/41 


°-3 



Hernandez ct al 
Hernandez et al 
Carpenter ct al. 
Hernandez ct al 



This paper 

Currie et al. (20081 
Smith et al. (20111 
Gorlova et al. (20071 




Notes: x Frequency of excess emission in this spectral range highly uncertain due to low number of stars. 
a No significant detections at 8^m in this spectral range. 



tics of iGehrelsl (|l986l ) as described for the 24/jm data in the 
preceeding section. 

These figures confirm that 8/im excess emission is rare 
amongst stars of ages ^5Myr. We do not see any evidence for 
a dependence of 8/xm excess frequency on the age of the clus- 
ter. As we only have upper limits for many of the clusters, 
particularly for the lower mass spectral bin, any dependence 
of 8/im frequency on age or spectral type is difficult to deter- 
mine. However we can state that the 8/mi excess frequency 
is always lower than the 24/^m frequency. This demonstrates 
that hot-warm dust is much less common than warm-cool 
debris disc emission. This in turn suggests that discs are 
likely to lie at a few to tens of au from the star, and not 
within 1 or a few au. This is consistent with the picture of 
debris discs arising from statis tical studies and studies of 
individual debris disc hosts fsee lWvattll2008l and references 
therein). 



6 CONCLUSIONS 

In this paper we have presented a study of the cluster NGC 
1960 (M36) using archival Spitzer MIPS and IRAC data. 
These data have been used to search for debris discs in the 
cluster. Our conclusions are as follows: 

(i) We have identified 38 targets as having significant 
24/mi excess. One target, bright star 14 (BD +34 1098), 
has very high levels of excess (i*24/i r phot ~85). Three tar- 
gets (including BD +34 1098) have significant excess emis- 
sion at near-infrared wavelengths. The SEDs of these targets 
have been examined and the emission has been determined 
to be consistent with second generation debris emission, as 
opposed to remnant primordial discs. 

(ii) The frequency and overall levels of excess emission 
observed in NGC 1960 are consistent with other clusters of 
similar ages (>30% of B0-F5 type stars, and >23% of F6- 
K9 stars have significant excess). The data from this cluster 
are consistent with a 'rise and fall' in debris disc emission 
amongst higher mass stars, and with a possible delayed peak 
in emission amongst lower mass stars. 

(iii) We find no evidence for a dependence of excess emis- 
sion on stellar multiplicity. We are however restricted in our 
examination of this relationship by a lack of completeness in 



the 24/im flux measurements for intermediate to lower mass 
stars. 



ACKNOWLEDGMENTS 

This work is based on archival data obtained with the 
Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology 
under a contract with NASA. This paper has made use of 
data obtained at the William Herschel Telescope, which is 
operated on the island of La Palma by the Isaac Newton 
Group in the Spanish Observatorio del Roque de los Mucha- 
chos of the Instituto de Astrofsica de Canarias. R.S. would 
like to acknowledge the support of the UK Science and Tech- 
nology Facilities Council (STFC). 



REFERENCES 

Abt H. A., 1983, ARAA, 21, 343 
Abt H. A., Levy S. G., 1976, ApJS, 30, 273 
Alonso A., Arribas S., Martinez-Roger C, 1996, A&A, 313, 
873 

Bryden G., et. al., 2006, ApJ, 636, 1098 
Canup R. M., 2004, ARAA, 42, 441 

Carpenter J. M., Mamajek E. E., Hillenbrand L. A., Meyer 

M. R., 2006, ApJL, 651, L49 
Carpenter J. M., Mamajek E. E., Hillenbrand L. A., Meyer 

M. R., 2009, ApJ, 705, 1646 
Chen C. H., Jura M., Gordon K. D., Blaylock M., 2005, 

ApJ, 623, 493 
Chen C. H., et al., 2005, ApJ, 634, 1372 
Chen C. H., Sheehan P., Watson D. M., Manoj P., Najita 

J. R., 2009, ApJ, 701, 1367 
Cieza L. A., et al., 2009, ApJL, 696, L84 
Cox A. N., 2000, Allen's astrophysical quantities. (New 

York: AIP Press; Springer) 
Currie T., Kenyon S. J., Balog Z., Rieke G., Bragg A., 

Bromley B., 2008, ApJ, 672, 558 
Currie T., Lada C. J., Plavchan P., Robitaille T. P., Irwin 

J., Kenyon S. J., 2009, ApJ, 698, 1 
Currie T., Plavchan P., Kenyon S. J., 2008, ApJ, 688, 597 
Dominik C, Decin G., 2003, ApJ, 598, 626 
Duchene G., 2010, ApJL, 709, L114 



16 R. Smith and R. D. Jeffries 



Duquennoy A., Mayor M., 1991, A&A, 248, 485 
Engelbracht C. W., et al, 2007, PASP, 119, 994 
Fazio G. G., et al., 2004, ApJS, 154, 10 
Feigelson E. D., Nelson P. I., 1985, ApJ, 293, 192 
Flaherty K. M., Pipher J. L., Megeath S. T., Winston 

E. M., Gutermuth R. A., Muzerolle J., Allen L. E., Fazio 

G. G., 2007, ApJ, 663, 1069 
Gehrels N., 1986, ApJ, 303, 336 

Gorlova N., Rieke G. H., Muzerolle J., Stauffer J. R., Siegler 
N., Young E. T., Stansberry J. H., 2006, ApJ, 649, 1028 

Gorlova N., Balog Z., Rieke G. H., Muzerolle J., Su 
K. Y. L., Ivanov V. D., Young E. T., 2007, ApJ, 670, 
516 

Grigorieva A., Artymowicz P., Thebault P., 2007, A&A, 
461, 537 

Hartmann L., Calvet N., Gullbring E., D'Alessio P., 1998, 
ApJ, 495, 385 

Hasan P., 2005, Bulletin of the Astronomical Society of 
India, 33, 151 

Hernandez J., Calvet N., Hartmann L., Briceno C, Sicilia- 
Aguilar A., Berlind P., 2005, AJ, 129, 856 

Hernandez J., Briceno C, Calvet N., Hartmann L., Muze- 
rolle J., Quintero A., 2006, ApJ, 652, 472 

Hernandez J., et al., 2007, ApJ, 662, 1067 

Hernandez J., Calvet N., Hartmann L., Muzerolle J., 
Gutermuth R., Stauffer J., 2009, ApJ, 707, 705 

Hernandez J., Morales-Calderon M., Calvet N., Hartmann 
L., Muzerolle J., Gutermuth R., Luhman K. L., Stauffer 
J., 2010, ApJ, 722, 1226 

Hillenbrand L. A., 2008, Physica Scripta Volume T, 130, 
014024 

Jeffries R. D., 2006. In Chemical Abundances and Mixing in 
Stars in the Milky Way and its Satellites, ed. S. Randich, 
L. Pasquini pp 163-170. Berlin: Springer- Verlag 
Jeffries R. D., Jackson R. J., James D. J., Cargile P. A., 

2009, MNRAS, 400, 317 
Kenyon S. J., Bromley B. C, 2002, ApJL, 577, L35 
Kenyon S. J., Bromley B. C, 2005, AJ, 130, 269 
Kenyon S. J., Bromley B. C, 2006, AJ, 131, 1837 
Kenyon S. J., Hartmann L., 1995, ApJS, 101, 117 
Klahr H., 2008, Nature, 52, 78 
Kroupa P., 2001, MNRAS, 322, 231 

Lagrange A.-M., Backman D. E., Artymowicz P., 2000, in 
V. M., P. B. A., S. R. S., eds, Protostars and Planets 
IV Planetary Material around Main-Sequence Stars. (Tus- 
can: University of Arizona Press), pp 639-672 

LaValley M., Isobe T., Feigelson E. D., 1990, in Bulletin 
of the American Astronomical Society Vol. 22 of Bulletin 
of the American Astronomical Society, Software Report: 
ASURV, The Pennsylvania State University, pp 917 — h 

Lefevre L., Marchenko S. V., Moffat A. F. J., Acker A., 
2009, A&A, 507, 1141 

Magnier E. A., Waters L. B. F. M., Groot P. J., van den 
Ancker M. E., Kuan Y.-J., Martm E. L., 1999, A&A, 346, 
441 

Makovoz D., Marleau F. R., 2005, PASP, 117, 1113 
Mayne N. J., Naylor T., 2008, MNRAS, 386, 261 
Moerchen M. M., Telesco C. M., Packham C, Kehoe 

T. J. J., 2007, ApJL, 655, L109 
Muzerolle J., Allen L. E., Megeath S. T., Hernandez J., 

Gutermuth R. A., 2010, ApJ, 708, 1107 
Naylor T., Jeffries R. D., 2006, MNRAS, 373, 1251 



Plavchan P., Werner M. W., Chen C. H., Stapelfeldt K. R., 

Su K. Y. L., Stauffer J. R., Song I., 2009, ApJ, 698, 1068 
Reach W. T., et al., 2005, PASP, 117, 978 
Rebull L. M., et al., 2008, ApJ, 681, 1484 
Rieke G. H., Lebofsky M. J., 1985, ApJ, 288, 618 
Rieke G. H., et al., 2005, ApJ, 620, 1010 
Rieke G. H., et al., 2004, ApJS, 154, 25 
Sanner J., Altmann M., Brunzendorf J., Geffert M., 2000, 

A&A, 357, 471 
Siegler N., Muzerolle J., Young E. T., Rieke G. H., Mama- 

jek E. E., Trilling D. E., Gorlova N., Su K. Y. L., 2007, 

ApJ, 654, 580 
Siess L., Dufour E., Forestini M., 2000, A&A, 358, 593 
Skrutskie M. F., et al., 2006, AJ, 131, 1163 
Smith R., Churcher L. J., Wyatt M. C, Moerchen M. M., 

Telesco C. M., 2009, A&A, 493, 299 
Smith R., Jeffries R. D., Oliveira J. M., 2011, MNRAS, 

411, 2186 

Stauffer J. R., et al., 2010, ApJ, 719, 1859 

Su K. Y. L., et al., 2006, ApJ, 653, 675 

Trilling D. E., et al., 2007, ApJ, 658, 1289 

Weidenschilling S. J., 1977, MNRAS, 180, 57 

Weidenschilling S. J., Cuzzi J. N., 1993, in E. H. Levy & 
J. I. Lunine ed., Protostars and Planets III Formation of 
planetesimals in the solar nebula, pp 1031-1060 

Wetherill G. W., Stewart G. R., 1993, Icarus, 106, 190 

Wyatt M. C, 2008, ARAA, 46, 339 

Wyatt M. C, Smith R., Greaves J. S., Beichman C. A., 
Bryden G., Lisse C. M., 2007, ApJ, 658, 569 

Yin Q., Jacobsen S. B., Yamashita K., Blichert-Toft J., 
Telouk P., Albarede F., 2002, Nature, 418, 949 



